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Abstract 

We review our previous work on the calculation of the vacuum expectation value of 
the stress-energy tensor for a scalar particle in the background metric of different types of 
spherical impulsive, spherical shock and plane impulsive gravitational waves. We get non- 
zero vacuum fluctuations only if we use the spherical impulsive wave as the background. 
We take a detour in de Sitter space for regularizing our result both in the ultraviolet and 
the infrared regions. 
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INTRODUCTION 

We know that the methods for the quantization of Einstein's theory of gravititation are 
not at their final stage at the moment. There are many models which incorporate this 
theory into the models for elementary particles. One can cite all the string, membrane, 
M models, etc. which are extensively studied I 1 . The final word is not still said in these 
models. 

In the absence of a generally accepted method for the quantization of gravitation, 
semi classical methods were used to extract information about the theory, which were 
very useful in the absence of a full quantization. We can calculate the fluctuations in the 
stress-energy tensor of a particle that propagates through universes described by different 
metrics. Extensive work in this direction was done in the seventies stressing the phenomena 
of particle production in these metrics . This work is described in the books written by 
Birrell and Davies / 2 , Fulling / 3 and Wald / 4 . Essentially we are confronted by a problem 
of a particle in an external field, as the metric contribution can be written as an external 
potential. The n-point functions where n exceeds two can be found by a simple application 
of the Wick's theorem in quantum field theory. Kuo and Ford show the calculation of the 
four-point function, thereby variance in the stress-energy tensor using this idea / 5 . 

In the eighties and nineties people applied the same method to problems in the pres- 
ence of gravitational waves, or to the effects of different topological structures in the 
space-time manifold. Vacuum fluctuations in the presence of cosmic strings is an example 
of such work / 6 . 

We were interested in the same problem in the presence of spherical waves created 
by the snapping of cosmic strings . We calculated the vacuum-expectation-value, VEV, 
of the stress-energy tensor, < T^ v >, for several types of spherical impulsive and shock 
waves, described by different warp functions. There is a general argument / 7 backed up 
by a explicit calculation / 8 which prohibits the polarization of the vacuum in the presence 
of plane waves. It was not known whether the same obstruction prevailed in the spherical 
wave case, though; so, the expression for < > was calculated / 9 for special cases of 
the Nutku-Penrose metric / 10 . We found that although a finite expression can not be 
extracted in the first order calculation, the second order calculation indeed gives a finite 
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result . Here we had to take a detour in the de Sitter space and take the appropriate limit 
to go back to the Minkowski space at the end of the calculation. 

The same method was applied / n to the spherical shock wave given by Nutku I 12 . In 
this case we were not able to extract a finite expression even in second order. 

We, then, applied our method to the impulsive plane waves. Relying on the Deser- 
Gibbons work / 7 ' 8 , we were not anticipating a finite result at this point. It seems that 
the method we proposed for this calculation passed this test. We could not get a finite 
contribution to the vacuum fluctuations if we had a plane impulsive wave propagating in 
the Minkowski space. If we have an impulsive plane wave in the de Sitter universe, though, 
we get a finite result which is proportional to the cosmo logical constant. We also find that 
the presence of impulsive pp waves does not change this result. 

In Section II we first calculate < T^ v > for an impulsive wave in first order pertur- 
bation theory for a special form of the warp function. In Section III we do the similar 
calculation for a specially simple form of the same function in second order perturbation 
theory. Section IV is devoted to the similar problem in the de Sitter background. Here we 
show how one can extract a finite expression for < T^ v > first by performing the calcula- 
tion in de Sitter space and then taking the appropriate limit to land in Minkowski space. 
In Section V we do the similar calculation for the shock wave and point to the absence of 
a finite term in the expression for < T^ v > even in second order in perturbation theory. In 
Section VI the similar calculation is repeated for an impulsive plane wave and it is shown 
how the proposed method fails to extract a finite expression for the vacuum expectation 
value of the stress-energy tensor in the background metric of a plane wave. This result 
is in agreement with the result given by Deser and Gibbons Z 7 ' 8 . We also show how this 
calculation can be extended to pp waves and give the result for this case. We end up with 
few concluding remarks. 

II. Spherical impulsive wave: first-order calculation 

The Nutku-Penrose metric / 10 for spherical gravitational waves formed as a result of 
snapped cosmic strings is given by 

ds 2 = 2dudv-u 2 \d( + vO(v)f(()d(\ 2 . 1 
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Here v is the retarded time, u is a Bondi-like luminosity distance, Q is the angle of stereo- 
graphic projection on the sphere. / is the Schwarzian derivative of an arbitrary holomor- 
phic function h{Q which describes the way the cosmic string is snapped and is called the 
warp function. Only one component of the Weyl tensor is not zero. It is given by 

* 4 = — /, 

u 

exhibiting the characteristic behaviour of an impulsive wave. We expect similar behaviour 
for the vacuum expectation value of the stress-energy tensor ,< T^ v >, for a massless scalar 
particle propagating in the background metric of such a wave if we get any finite value for 
this expression. 

We choose a particular form for the warp function h to perform an explicit calculation. 
Our experience with the different forms we chose / 9 shows that the form of the expression 
for < T Mi , > is not sensitive to the actual choice for h. The presence of the Dirac delta 
function and the presence of a finite expression depends only on the general properties of 
the metric. The functional form of the warp function only effects the dependence on £, (. 
This fact allows us to choose a special form for this function that will make our calculations 
as simple as possible. 

Our first choice corresponds to a "rotating" string and is given by 

h=(C) 1+ie 3 

which gives 

/(C) = ^ 4 

in first order in e. With this choice we rotate the shear of the metric by 7r/2 compared with 
the real exponent case. We take the rotating string because the calculations are somewhat 
simpler in this case. Here e is a small number which can be positive or negative, and will 
be used as the perturbative expansion parameter. 
First we go to real coordinates 



and write the metric, up to first order in e, for v > 0, 

ds2 = 2dudv ~ ( u2 ~ (frw) dx2 -{ u2 + (frw) d y 2 - jT y ^ {x2 - y2)dxdy - 
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The d'Alembertian operator □ times y/—g , written in the background of this metric reads 



n o d 2 n 8 2 2 8ev 
2-u 2 -— — + 2u— 



( f) 2 f) 2 \ „ „ B 2 1 

. 7 



.,„, d 2 d 2 \ 2 y2 d 2 
dx 2 dy 2 J dxdy 



dudv dv dx 2 dy 2 u(x 2 + y 2 ) 2 
For technical reasons we first couple a massive scalar field to this metric and take the 
limit mass going to zero at the end of the calculation. We use conformal coupling. The 
Ricci scalar is zero in the Minkowski case. The form of this coupling only matters in the 
de Sitter case. 

We first compute the Green's function by summing the eigenvalues of the operator 

G F = J2 t^MM^l 8 

A 

where 

L<j>\ = 9 

L = > / == ^[n + m 2 ]. 10 

We calculate < T^ v > by taking the appropriate derivatives of Gp. Since our computation 
is in first order in e, we expand L, <p\ and A and write 

(L + eLi)(0o + e(j>! + ...) = (A + eAi + ...)(0 O + e^i + •••)• H 

This gives us 

L (f) = A o 0o, 12 
Li0 o + ^o0i = Ai0 o + A o 0i, 13 

where 

/ d 2 d d 2 d 2 \ 

\ anew at> oar ay 2 / 
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e iRv e ik 1 x e ik 2 y e ^§ e i?pL 



?0 = 



Uy /\2R\W 



15 



X = k 2 + kj-K 
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8v 

Li — — 

u 



xy 



d 2 



(x 2 + y 2 ) 2 \dx 2 dy 2 J (x 2 + y 2 ) 2 dxdy 



2 2 

x z — ?r 
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Ai = (00,-^100) = 0. 



We write 



01 = 0o/- 

The form of the operator L\ given in eq. 17 suggests the ansatz 

/ = vf!(x,y,u) + f 2 (x,y,u) 
for /. This ansatz gives us two coupled equations 



L'h = 



8 



xy 



(x 2 + y 2 ) 2 
df 



u 

L'f 2 = -2u 2 ^ 
ou 



{-k 2 + k 2 ) + 



(x 2 -y 2 )k x k 2 



(x 2 + y 2 ) 2 
iK 2im 2 u 2 



1T + 



where 



R 

d\ d 2 



A, 



3 fa , a \ a 2 a 2 

L = 2u 2 iR— - 2i fcx— + k 2 — - ^ - j—. 

au \ ax ay j ax z ay z 
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Here we take mass only as an infrared regulator and from this point on equate it to zero 
everywhere except in the zeroth order solution 4>q. We will take the limit where it goes to 
zero at the end of our calculation. 

At this point we see that it is convenient to make the change of variables and use 
z = x + iy,z = x — iy, and s = ^. In terms of these variables we get 



L' = -2iR^- - 2i 

as 



(h + ik 2 )^- + (h - ik 2 )-^- 

dz dz 



a 2 



dzdz 



L'fr = -2is 



(k x +ik 2 ) 2 (h-ik 2 y 
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with the immediate solution 

h = R(lnz - Inz) - ( (fel+ ; fe2)s - 26 

h = - i) - Ulnz - Inz) + k( -* - ^ ~ Z 1 . 27 

If we go back to the real variables the solutions read 



A = 2Ritan~ 1 - - 2is^ ^ 28 

__l_ rt,A 



f 2 = -(—--)tan ~+ ^ + k 2 [2(hx + k 2 y)(tan - - y) 

+ (k t y - k 2 x) ((Inx 2 + y 2 ) - 2x)) 29 
To form the Green's function Gp, we have to calculate 

G f = r *, r * 2 r iK r ffl »wMM, 30 

</-oo J-oo J-oo J-oo K k 1 k 2 

where 0q and 0^ are the zeroth and the first order solutions respectively. We have to 

multiply fi and f 2 by the operator 

f°° f°° [°° [°° 1 f 00 dR -iK(ss') 

L = % I da I dk\ I dk 2 I dK- — ^- — — - / — r-=-e 2R 

Jo J-oo J-oo J-oo (87r 2 )(27r)Wi_ 0O 2| J R| 



xe 



iR( v - v ') e ik 1 (x-x') e ik 2 (y-y') e ic t (K-k 2 1 -k 2 2 ) e ^-(^-j r ) gj 



At this point we are looking for finite contributions to the < T Mi , >. To calculate them we 
have to first calculate the Green's function, and then differentiate it to obtain the vacuum 
expectation value of the stress- energy tensor. Since we expect a non zero contribution, if 
any, only from < T vv >, taking the coincidence limit in x and y before we calculate Gf 
will not change the end result because we have to differentiate with respect to v and v' 
only. If these terms are zero in the coincidence limit in x and y for the Green function, 
they will be also zero for < T^ v >. 

A straightforward calculation shows that if we take the coincidence limit for x and y, 
i.e. take x = x' , y = y' , we find 

(u — u ){v — v ) (u — u')[v — v'Y X 
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Here A, B are numerical constants. We could not extract a finite portion out of these 
expressions in the coincidence limit. 

One can show that if we go one order higher, we may get finite contributions if we 
take a detour in de Sitter space. Since this calculation is quite involved, we first use a 
simpler form of the warp function and demonstrate how this method works. 

III. Spherical impulsive wave: second-order calculation 

In this section we take the warp function h = e Q< > , which results in / = resulting in 

a metric 

ds 2 = 2dudv - - (dx 2 (2u - va 2 Q{v)) 2 + dy 2 (2 + va 2 Q{v)) 2 ) . 33 
The d'Alembertian operator in this metric reads 

□ = -£=Wy/=g)d v = 2d^d u + ( — L- + — i-j- ] d v 



-9 " \u-^ u+Y 
+ ^! [ __L L_ ] du l -^d 2 x 1 -^d 2 v . 34 
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for the exact operator. We multiply this expression , eq. 34, by ^f^g which is equal to u 
and expand the operator up to second order in a 2 : 

L H = 2u 2 d u d v + 2ud v -dl-d 2 y -— (d 2 x -d 2 )- a 4 ( V -d u -—^ + §4 ( d * + d2 y)) • 35 

Here, again, we will first compute the Green's function for this operator, and then 
differentiate it to get the vacuum expectation value of the stress-energy tensor. We obtain 
the same zeroth order solution as given in eq. 15 and make the same ansatz (p^ = </>(°) / 
since the first order operator is of the same form as previous case. We see that / just 
modulates the zeroth-order solution, and does not essentially change it in a radical manner. 

/ obeys the differential equation 

L 2 f=-(k 2 -k 2 ) 36 
u 

where L 2 is defined as 

I - > -2iR— -2i(k — + k — \ -— - — -2-^— V— — } 37 

' ds \ 1 dx 2 dy J dx 2 dy 2 dsdv R dv J 



Here the modes fci, K 2 , -R, K are as given in eq. 16, s = K We make the ansatz for / 
similar to the one given in eq. 20. / = vfi(s,x,y) + f2(s,x,y). This yields equations 
similar to the ones given in eq.s 21 and 22. We get 

L'h = - 38 
L'h = (^ + '^)h 39 



where V is as defined in eq. 23. These equations are simply integrated over with the result 

To get the Green's function at this order we operate on this function by L given in eq. 31. 
The result is seen to yield 

where 

[ ] = («-« / )(^-^)-^((^-^) 2 + (2/-2/ / ) 2 ), 41 

and Ai,A 2 , A 3 are constants. We see that this result is of the Hadamard form. We find 
that all these terms have the same type of ultraviolet singularity as the flat part. We could 
not find a finite part of this expression. 

If we go one order higher, we end up with the differential equation 

L 2 g = v 2 (^ + 3s 2 (k 2 + k 2 ) + (k 2 - k 2 ) 2 

f-s iKs 2 (k 2 -k 2 ) 2 s 3 iK(K 2 -k 2 ) 2 s 4 \ 
+V {^ + ^R- 4R 2 + 24R3 ) 42 

when we make the ansatz 4> 2 = 4>og. We take 

g = v 2 g 1 (x, y, s) + vg 2 (x, y, s) + g 3 (x, y, s). 43 
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Going through similar steps we find 

.2 A<S(1.2 



9i 



-s- is\k 2 + k 2 ) | {k 2 -k 2 ) 2 s A 
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92 = 



-i3s 2 13Ks s is 



2/? 32i? 2 

4 
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8R 



+ 



12R 2 32R S 



((k 2 -k 2 ) 2 + K(k 2 + k 2 )) 



160.R 4 



(k 1 — k 2 ) , 
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3s 2 i55Ks 3 



93 = 



8R 2 16R 3 

is 5 K 



^{iki-k 2 r + K(ki+ki) + 1 ^f 



(3(kl-kl) 2 + K(kl + k 2 2 ))- 



K 2 s 6 (k 2 -k 2 ) 2 
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32(1R 5 v v 1 2/ \ l > zjj 1920.R 6 
We see that when the Green's function is calculated using these functions we get no 
infrared divergences for g\. When we use g2, we have to perform integrals of the type 



Jo 



da 



a 



-exp 



-i(u — u')(v — v') 



uu a 
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which results in a soft infrared divergence. If we use an infrared regulator mass m 2 , we 
get terms that go as log ((m 2 (u — u'){v — v')). m 2 cancels out when derivatives are taken 



to obtain < T^ u >. 



The infrared divergences are more severe for (73. There we have to perform integrals 
of the type 



daexp 



—1- 



(u — u')(v — v') 



uu a 
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which are linearly divergent at the upper limit for the massless case. If we use a massive 
field as an infrared cut-off, we get terms that go as ^ as m tends to zero. This term 
multiplies the whole expression and does not drop out on differentiation. 

IV. Going to de Sitter space 

Here we treat the same problem for an impulsive wave in de Sitter space. It is known 
that an exact solution can be found for an impulsive wave propagating in de Sitter space 
if we multiply the Minkowski space solution by the conformal factor (l + ^p) Z 13 - This 
is reflected to the Gp by factors multiplying the expression from both sides. 



G% = 1 + 



Auv\ 



GM(x,x>)[l + 



Au'v' 
6 
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where Gp and Gp are the de Sitter and Minkowski space Green functions respectively. 
Note that 

1 Auv\ / AuV\ / AUV\ 2 A. 



12. 

where U = , U = Using this expression we see that for Gp^ 1 ' we get terms 

in the coincidence limit that give a well-defined finite expression for < T$ >. Here the 
regularization is performed using the point-splitting method, exactly in the same way one 
would do for the zeroth part of the de Sitter universe. One sees that for the warp function 
we studied in Section II, we get a finite contribution. If we use h = (C) 1+s+te we get 

< T - >= ( Sl ° 3{x2 +y2) ~ 4ftan_1 ! ) ■ 61 

This result is finite, contrary to the case for propagating in Minkowski space. This is no 
contradiction to our previous results, since this expression vanishes as A goes to zero. One 
can also show that ' 1A for h = ( J we get 



< 



All of these calculations are in first order perturbation theory. We saw in Section 
III that there are problems in the second order calculation due to the appearance of 
infrared divergences for the special warp function used. One can show that these infrared 
divergences are generic for the Nutku-Penrose impulsive metric Z 15 in second order and all 
can be tamed to give a finite expression for < T^ v > in the following way. 

The solution for the same type of wave in de Sitter space has the multiplicative factor 
^ multiplying a finite expression. This fact suggests taking the curvature of the de Sitter 
space A proportional to the infrared parameter m 2 . These two terms both have the same 
dimensions, as were introduced only as a technical aid to get rid of the ultraviolet and 
infrared divergences. Our end result should be independent of them. This is achieved if we 
take them to be proportional to one another, and let them go to zero at the same rate. The 
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finite proportionality constant that may appear is absorbed in the perturbation expansion 
parameter, a 2 for the case studied in Section III. 

Studying this particular case, one sees that applying the operator L as given in eq. 
31 on the expression (73, eq.(46), results in divergent terms of the form 

(s — s') 2 m 2 (s — s') 2 m 2 

in Minkowski space. If we go to de Sitter space, we have to multiply this solution by the 
factor ^(u — u')(v — v'). This operation gives rise to a finite expression when both A and 
m 2 are taken to zero at the same rate. 

The result of this exercise is that by performing the calculation in de Sitter space and 
taking the limit where the curvature of this space goes to zero at the same rate as the 
infrared parameter, we indeed get 

6(v) 

< T vv >oc -V 54 

which is proportional to the sole non-zero component of the Riemann tensor , as anticipated 
in the first place. 

In the next section we see that this property is very special to the impulsive spherical 
waves and is not shared with the spherical shock or impulsive plane and impulsive pp 

waves. 

V. Spherical shock wave 

In first order perturbation theory both the impulsive and the shock wave cases showed 
similar behaviour. The two solutions are essentially different, though. The shock wave 
solution has a dimensional constant which is lacking in the impulsive wave solution. Since 
in quantum field theory , models with dimensional and dimensionless constants belong 
to different classes, we thought similar distinction between these two models may exist. 
Relying on these motivations we planned to check whether there is a qualitative difference 
between the two solutions exhibited by their behaviour at higher orders . 

We will show that the infrared divergences which may be cancelled via a detour in de 
Sitter space are absent in the shock wave calculation. Another point of difference is the 
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importance attributed to the homogenous solutions in these two cases. The homogenous 
solutions give just the free Greens function for the impulsive case, whereas they result in a 
totally different contribution to the Greens function in the shock wave. This is an artifact 
of the presence of a dimensional coupling constant in the latter case. What may be more 
interesting is the fact that just the contribution of the first order calculation contributes 
to < > in de Sitter space. The higher order terms cancel out when the VEV of the 
stress-energy tensor is computed. 
We start with the metric Z 12 , 

ds 2 = 2Pdudv + 2uP c d(dv + 2uP^d(dv - 2u 2 d(d(. 55 

Here we use the same variables as those given for the impulsive wave; i.e., u is a Bondi-type 
luminosity distance,?; is a timelike coordinate that can be regarded as retarded time and 
C and ( replace the angles on the stereographic projection of the sphere. P is defined by 
P = y^-jj, where h is an arbitrary function of the argument £ + gv@(v). g is the dimensional 
coupling constant, with dimensions of mass and is the Heavyside unit step function. 
Nutku Z 12 shows that the identifying factor in this metric is one component of the Weyl 
tensor, 

*4 = -^p[PP<C<-P«P<;]9(v)- 56 
When we write the d'Alembertian operator in the this background metric, we get 

\/~~gU L 2u q u q v + u q v + u CC g u u p Q u 2 U P du 

( d d \ d d 2 

+2u[P 7 —+P c -=) — - 2P—= 57 

In our particular case we take h = ((+gv Q(v )) 1+J<5 , where 5 « 1 and is the expansion 
parameter. 

We expand the operator L in powers of 5. 

L = L + 5Li + 5 2 L 2 + ... 58 

If we use real variables, ( = ^7=^, we get 

« 2 d 2 n d d 2 d 2 
ouov ov ox z oy z 
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u = 



2u 



d d \ d _ lV ( d 2 d 2 

Vt{ x— I — tan — I t^-tt + 



x 2 _|_ y2 y Q x Qy J Q u 



x \dx 2 dy 2 



60 



L 2 = 



1 (d_ 

x 2 _|_ y2 ^ ^ 



2 <9 2 _ _i y / 9 <9 \ <9 

+ 2tcm -h/a 

x \ ox oy J ou 



+ 



d 2 d 2 
+ 



dx 2 dy 2 J ' 

We expand both the solutions and the eigenvalue in terms of <5, 

= 0o + <tyi + 5 2 2 + A = A + <5Ai + S 2 X 2 + .... 
Up to second order we get the coupled set of equations 

L (p = A o 0o, 
Li0o + Locp! = Ai0 o + A o 0i 
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L 2 (f>o + £i0i + -^002 = A 2 0o + Ai0i + A o 02 65 
We make the ansatz 0i = f4>0i 02 = ^0o- The ansatz for 0i results in the equation 

a 2 a 2 



a 2 a (i^l\_2i(k— k — 

dsdv ds \ R dv ) \ 1 dz 2 dy J dz 2 dy 2 



f 



2i 



(hy-k 2 z) - i) - (fc? + A*)f an" 1 (|) 



66 



^2 _j_ ^2 

We recognize that when we plug these ansatze into our equations, we get inhomogenous 
equations where on the right hand side the variable v always appears in the combination 
z = x + gv for v > 0. v is not an independent variable in / and h. Thus we substitute 
to new set of variables u,z,y. This necessitates replacing all v derivatives by g-^. This 
makes an expansion in powers of g possible. We take 

/ = fo(z,y,u) + gfi(z,y,u), 

h = h (z,y,u) + ghi(z,y,u). 67 

0o is given by, as in the previous cases 

exp[i(fcia: + k 2 y + Rv - ^)] 

<f>0 = 7= — • Do 



u^/\R\(2tt) 2 
14 



where K,k±,k 2 , R are the seperation constants which act as eigenfrequencies to be inte- 
grated over to find the Greens Function. For z = x + gQ(v)v, s = ^ we find 



£ /o = I 
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where 



£1/0 + £i/o = 



C = -ZiRjr- -2t(k 1 — + ki- 
ds \ az ay 



£ -2-^— — — 

1 dsdz R dz 



«9 2 



Qy2 Q z 2 



I = 



2i 



(z 2 + y 2 ) 



(hy - k 2 z) 



iKs 
~2R 



1 



{k{ + k^tan' 1 - 
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We see that if we use the variables z + iy and z — iy the right hand side can be written as 
a sum of functions of this two variables. We also write the differential operator on the left 
hand side in terms of these new set of variables and after some algebra we find 



fo = 



Ks i 
JE + 2 



In 



z + iy 
z-iy 



K + kf + k^ 



1 



(&4 + ik 2 ) 



[{z + iy)ln(z + iy) - (z + iy)] 



~Ti — ~~it~\ tt z " ^) /n <> - w) - - iy)] ) 

(ki-ik 2 ) J 

I iK K 2 s\ ( ln(z + iy) ln(z — iy) 
h = \2R + 8R 2 ) V h + ik 2 kl 
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ik-2 



K 



+ 



K 



((&i - ik 2 ) 2 ((z + iy)ln(z + iy) - (z + iy)) 



AR \2(k 2 + k 2 ) (kj + k 2 ) 2 , 

-{ki + ik 2 ) 2 ((z - iy)ln(z - iy) - (z - iy))) . 
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The Green function in first order for the warp function chosen above is given in refer- 
ence 16. The expression given in this paper was calculated by a different method and also 
includes the contribution coming from the homogenous solution. Here we calculate them 
seperately . We can show that we do not get a finite contribution to the vacuum expecta- 
tion value of the stress-energy tensor. When the contribution only from the inhomogenous 
equations (eq.s 69,70) is taken, the Green function is given by a product of two factors one 
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is which always equals zero when the limit in x = x' , y = y' is taken without equating u, 
v to u',v' respectively. 

In second order in 5, we can reduce the differential equation to the system 



C h = M , 76 
Coh! +dho = Mi, 77 



where 



d f d d \ d 2 d 2 

£ = -2iR—-2i(k 1 — + k 2 —)- — -—^, 78 
os \ oz ay J oz z ay 1 

C - 2 ° 2 h — — 79 

1 dsdz R dz 

Here s = ^. M and Mi are given as 

-2K ^ l °9( z2 + y 2 ) _ ^ {kiy _ hz) + tan -ii ihz + hy) 

Xiil -^r) (z 2 +y 2 )K -2 tan ~z 



z 2 + y 2 



_ 2iyK f %Ks\ _ x y 2k x ( iKs K 2 s 2 

1 ~ (z 2 + y 2 )R V ^Rj tan ~z~ (z 2 +y 2 )R V ~ ~8W 

- (Zkxtan- 1 ^ - k 2 log(z 2 + y 2 )^j 

ff-iK 2 s 5K\ _ x y i(hy-k 2 z) f iKs K 2 s 2 \\ 
{\jR 2 ~ + 4Rj tan ~z ~ (z 2 + y 2 )R \ ~R~ " W ) ) 

+ {{k\ - k 2 2 ) (2ztan- lV - + ylog(z 2 + y 2 ) - 2y) 

-2k x k 2 [zlog^z 2 + y 2 ) - 2ytan~ 1 ^ - 2z^j 

( 3Ki ,y 3 ., , iKs.\ 3i / iKs* 

x tan- 1 ?- - —————[k xy - k 2 z)(l - —) - 1 - 
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In these expressions we took the 'mass-shell' condition, which is imposed in the calculation 
of the Greens function; i.e. we set k 2 + k 2 equal to K. One can check that after we perform 
the K and fci, ki integrations the effect of these two expressions are exactly the same. 

We see that, contrary to the impulsive wave calculation, in both M and Mi, there are 
no terms that are independent of z and y except a single term which is proportional to K. 
To be able to obtain terms in < > that diverge as the infrared parameter goes to zero, 
we need inverse powers of R which are not multiplied by K or k 2 , k\. Each inverse power 
of R means a higher infrared divergence, resulting in terms that go m 2 , 1, logm 2 , ^7, 
etc. ..as the power of is increased, whereas each power of K, k\, k 2 means one lower 
order in the same divergence. In free space the power of m is zero. There is no infrared 
divergence. 

In reference 15, we generated these divergences at second order and then cancelled 
them with the cosmological constant of the de Sitter solution. Our mechanism for obtaining 
these infrared divergences was as follows. We isolated the —2iRJ^ in the operator £ horn 
the others and equated it to the term which did not contain z or y. 



where c can be a function of v but not that of z and y. Then h = has an extra power 
of 1/R compared to the other terms. The second iteration gives us h[ oc s 3 /R 2 . Such 
a term will induce 1/m 2 factor in the expression for the Greens Function , Gf, and this 
infrared mass will be retained in < >. 

For the shock wave solution all the terms in M and Mi are either functions of z and 
y or are multiplied by K. We can not isolate a part of the operator jCq and equate it to a 
single term on the RHS. We may still study the shape of solutions. Let us make the ansatz 



—2iR—h' = cs 
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h = s 2 h 00 + sh 01 + h 02 . 
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We find that 
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hoi = — rrr- f tan - ) 85 



etc. These solutions show that we are not getting powers of R in the denominator with no 
powers of K in the numerator. This fact makes it clear that there is no way of acquiring 
in the Green function, since we need an extra R~ 2 factor, which necessitates at least one 

TV" 2 TV" TV" 2 

power in the first integration. On RHS only the combination and exist . gives 
exactly the singularity structure as the free case, and jj? gives a logarithmic divergence 
which is cancelled in the < > calculation. 

At this point we note that we can find solutions of equations 76 and 77 even if M 
and Mi are set to zero. These are the homogenous solutions of the problem which give a 
non trivial contribution for the shock wave calculation. These solutions of the homogenous 
equations generated the only non-zero part of the Green function when we took x = x',y = 
y' in the first-order calculation. Since Mq and Mi are independent of v we can assume a 
powers series expansion in v for a chosen order in g. For the sake of illustration we take a 
solution in third order in g and write the expansion as 

^' 3) = 9^f^ + - 2 /& 3) + vf^ + /& 3) ). 86 

Here f^jj (s, z, y) has dimension zero, f^lf^i 8 , z i y) nas dimension minus one, etc.. Inverse 
powers of v are excluded by the regularity at v = 0. One can show that taking powers of v 
higher than that of g do not give results that differ from the free case. A similar expansion 
in the impulsive case would go as 

f£' 3) = (v*R*fi]? ) +v>R 2 f^ + ...) 87 

(1 3) 

where etc., have the same dimensions as above, since the only free dimensional 

parameters are v and R. This gives the result of the flat , Minkowski, space Green function. 

Keeping track of powers of v we get a system of four equations. We note that the first 
of these equations 

£ /& 3) = 88 



has a solution for any function 



F = F(^-(ki±ik 2 )-(z±iy)). 89 
18 



We can also show that the singularity behaviour of the Greens Function is independent of 
the form of F. At the end we get, for the worse infrared poles the expressions 

H V (u - u')(v - v') J 

= 2- 2 ( ^gW±gg^ log(2^ (ll - „> - V))) , 91 

\ (it — w') 4 m 2 / 

H V (w - u ) y 

Here c; are functions of a; and y, depending on the form for F used, m is the infrared 
mass. If we use a linear functionfor F, is proportional to y or x. 

Upon symmetric differentiation with respect to v and v', the terms with ^jy and ^ 
vanish. We find a finite contribution only if we go to the de Sitter space, i.e. multiply by 
the factor (1 + A f £ )(l + In this case we get 

< T vv >= g 3 (hAvQ(v) + f 2 A 2 uv 2 e(v)) 94 

which goes to zero with A when we go back to Minkowski background. In this expression 
fi are regular functions of x, y. One can also show that any terms with less diverging 
powers of m~ 2 in G^' H ^ do not give a finite contribution even in de Sitter background. 

Here we tried to show that two qualitative differences exist between the shock and the 
impulsive wave solutions proposed by the same group / 10 ' 12 . In the shock wave solution 
the infrared divergences which may be used to tame the ultraviolet divergences to result 
in finite contributions to < T vv > are absent in second order perturbation theory. We 
can not find finite contributions to < T vv > in Minkowski space. If we go to de Sitter 
space, though, we get a finite contribution which is proportional to function, which is 
the signature of a shock wave solution. 

The homogenous solutions, in the shock wave, give contributions to the Greens func- 
tion expression which are different from the free case. These solutions also give a finite 
contribution to < T vv > in de Sitter space. The presence of these nontrivial solutions is 
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only due to the dimensional coupling constant. The presence of g in the expansion makes 
it necessary to have an extra power of in the solution which results in a nontrivial term 
in Gf- 

VI. Plane impulsive wave 

As the last example here we take the metric describing an impulsive plane wave Z 17 , 

ds 2 = 2dudv - \d( + q cc vO(v)d(\ 2 . 95 

If we take q = g^- we get a plane wave. If the power is higher than quadratic we get pp 
waves Z 18 . The d'Alembertian operator in this metric is written as 

If we expand up to second order in the coupling constant g, we get 

L » 2d u d v - 2vg 2 d u - (1 + 3(vg) 2 )(d 2 x + d 2 + 2vg(d 2 x - d 2 y )). 97 

The zeroth-order solution is similar to the previous cases, eq. (15) resulting in a Green 
function that goes as 

A 98 



(u — u')(v — v') — \[{x — x') 2 + (y — y') 2 ] 

for constant A. We expand the solution in powers of g and take the first order solution as 
cpi 1 ) = f(p . It is straightforward to solve for / and we get 



f (M - fc 2> 

1 2iR 



i Ku 
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For the second order solution we again take 4>^ = (p^h. Here h = v 2 hi{x,y,u) + 
vh2(x,y,u) + hs(x,y,u). A straightforward calculation gives us 

hi = ^(k 2 1+ k 2 )u-^(k 2 -k 2 ) 2 , 100 

= -^(f - m + k 2 )) - ^ m - kif + m 2 + k 2 )) + K[kl -* 2?u \ ioi 
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hs = - 3(fc? + fc 2 2 )) + (3(^ - fc 2 2 ) 2 + 3K(fc 2 + k 2 ) - K 2 ) 

Here we see that a peculiar thing happens. We almost get a finite expression, but we have 
one power of (u — u') too many in the denominator. If we go to de Sitter space to cancel 
both the ultra-violet and infrared divergences, we see that the finite part of < T vv > goes 
as 

A 2 

< T vv >oc -2 — Q(v) 103 

which is finite only in de Sitter space. One power of the curvature cancels with the infrared 
parameter since we take A oc m 2 , but the remaining power takes the contribution to zero 
when we go back to Minkowski space. 

This result which is in accord with general arguments of Deser Z 6 ' 7 is a check that our 
method does not contradictany known results. 

One can show that this result does not change in the presence of a pp-wave background. 
Whether a wave is plane or pp type depends only on the form of the function q(Q in the 
metric. The general behaviour of the expression for the vacuum expectation value of the 
stress-energy tensor does not depend on the form of the function q. This form only changes 
an overall factor which can not decide whether the whole expression is finite or null. The 
same behaviour was already seen in the different warp functions we have used for the 
spherical wave. 

CONCLUSION 

Here we reviewed the work on the vacuum fluctuations for the stress-energy tensor for a 
scalar particle in the background of different impulsive and shock wave metrics. We found 
a finite fluctuation only for the spherical impulsive wave metric, whereas our result was 
null for the spherical shock and plane impulsive waves. We attribute the difference in the 
two cases to the presence of the dimensionless and dimensional constants. 

If we calculate the fluctuations for a conformal metric, fluctuations should be absent I 2 . 
We first perform perturbation theory about the Minkowski space, and our perturbations 
are not strong enough to overcome the restrictions imposed by conformal symmetry. If we 
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go to de Sitter space, and perform perturbation around that metric, we do not have this 
obstruction. We always find finite fluctuations in that metric. Note that de Sitter space 
also tames our ultra-violet divergences. It turns out that if we have dimensional coupling 
constants, we have more severe ultra-violet divergences which are tamed only with having 
higher powers of the curvature scalar of de Sitter space, multiplying our expressions for 
the fluctuations. This happens in the two latter metrics, plane and pp waves, we have 
studied . The resulting infrared divergences are not strong enough to cancel these scalar 
curvature factors to give us a finite result. 

In the spherical shock wave solutions, we could not generate such divergences in the 
inhomogenous case. For the homogenous solutions, the ultra-violet divergences could not 
be tamed by going to de Sitter space . Cancellation of the infrared divergence , say in 
eq.92, by cancelling it by the A term was not sufficient to obtain a finite expression. 

The presence of non-trivial homogenous solutions was another novel feature of the 
shock waves. It was amusing that only one of them survived in the expression for < T^ v >. 
If many of them survived, we would have an ambiguity in defining this quantity. It was 
also interesting that no matter how many orders in the coupling constant one expands, 
one gets one and only one surviving term, the one obtained already in first order. 
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